Fiber bundle model with stick-slip dynamics 
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We propose a generic model to describe the mechanical response and failure of systems which 
undergo a series of stick-slip events when subjected to an external load. We model the system as 
a bundle of fibers, where single fibers can gradually increase their relaxed length with a stick-slip 
mechanism activated by the increasing load. We determine the constitutive equation of the system 
and show by analytical calculations that on the macro-scale a plastic response emerges followed by 
a hardening or softening regime. Releasing the load, an irreversible permanent deformation occurs 
which depends on the properties of sliding events. For quenched and annealed disorder of the failure 
thresholds the same qualitative behavior is found, however, in the annealed case the plastic regime 
is more pronounced. 

PACS numbers: 46.50.-|-a,62.20.Mk,81.40.Np 



There is a large variety of systems which undergo con- 
formational changes when subjected to external mechan- 
ical loads. The dynamics of spatial rearrangements can 
usually be described by the stick-slip mechanism, i.e. 
when the local load exceeds some threshold value, sub- 
units of the system increase their length leading to re- 
laxation. The stored length which can be activated by 
external loading, typically arises due to the existence of 
frictional contacts or chemical bonding between subunits. 
Several examples of such systems can be mentioned from 
the molecular scale organization of spider silk [l| , through 
the chains of magnetic beads in magneto-rheological flu- 
ids on the meso-scale 0, [1] , to the wire nets used to pro- 
tect roads from rockfalls in mountains 0] on the macro- 
scale. Experiments have revealed that the surprisingly 
high fracture toughness of spider silk is partly caused 
by the presence of blobs of protein molecules with a 
folded hairpin structure, which get unfolded under exter- 
nal loading [l[. A similar mechanical response has also 
been observed for biological tissues composed of intercon- 
nected bundles of fibrils 0| . In magneto-rheological fluids 
particles of permanent magnetic moment aggregate and 
form chains aligned with the external field, which then 
modify the rheological properties of the fluid. Stretch- 
ing and bending chains of particles of micrometer size by 
optical tweezers revealed that the anisotropic magnetic 
interaction and the frictional contact of particles result 
in a noisy response, i.e. particles undergo subsequent re- 
arrangements increasing the length of the chain and re- 
ducing the reactant force 0, Q • A similar mechanism is 
exploited on the macro-scale in wire nets, which cover 
steep walls in mountains protecting roads from rockfalls. 
In order to dissipate the kinetic energy of falling boulders, 
the wires form rings which can slide and collapse without 
breaking resulting in a large dissipation but keeping the 
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integrity of the net [J| . 

In the present paper we introduce a micro-mechanical 
model which captures the main ingredients of the me- 
chanical response of system which undergo conforma- 
tional changes with stick-slip mechanism. Our model 
construction is based on fiber bundle models 0, 0, d, 0, 
lol . m extended in such a way that fibers undergo slid- 
ing events which gradually increase their relaxed length 
when the local load exceeds some threshold values. Af- 
ter a large number of sliding events, fibers may also fail 
under a large enough external load. We derive the con- 
stitutive equation of the system and show analytically 
that the sliding mechanism leads to macroscopic plas- 
ticity of the bundle and permanent deformation remains 
after the load has been released. We explore the case of 
both quenched and annealed disorder of sliding thresh- 
olds. 

Our model consists of iV parallel fibers which have 
identical clastic properties characterized by the Young 
modulus E 

@, Si, EE, El- Under an increasing ex- 
ternal load (T parallel to the fibers' direction, the fibers 
exhibit a linearly elastic behavior until the local defor- 
mation Si reaches a threshold value ej,j, i = 1,.. .,N. 
The novel element of the model is that at the threshold 
the fiber does not break, instead it suffers sliding, i.e. its 
relaxed length increases until the fiber becomes capable 
to sustain the remaining load. The sliding threshold of 
fibers Eth is a random variable with a probability density 
pisth) and a distribution function P{eth) defined over the 
domain [e™",ej^°^]. The sliding event is instantaneous, 
i.e. it docs not take time. After the relaxation the fiber 
can be loaded again which may lead to further slidings. 
The total number of allowed sliding events (steps of ex- 
tension) kraax IS an important parameter of the model. 
When the fibers can slide more than ones kmax > 1, 
they either keep the initial sliding threshold for all slid- 
ings (quenched disorder), or can obtain a new threshold 
value each time from the same probability distribution 
(annealed disorder). The constitutive behavior of single 
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FIG. 1: Response of a single fiber for quenched (a) and an- 
nealed (b) disorder. When the failure threshold is reached 
the load on the fiber drops down to zero. The threshold can 
be constant (a), or new ones can be drawn from the same 
distribution (fe). 
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FIG. 2: {Color online) Constitutive behavior of the bundle 
with exponentially distributed quenched failure thresholds. 
a{e) tends to an asymptotic linear response preceded by a 
more and more horizontal plateau when kmax increases. The 
intersection of the asymptotic straight lines with the horizon- 
tal axis indicates the value of ej""^ obtained when unloading 
the system. For annealed disorder the cases of kmax ~ 5 (red, 
dotted line) and 10 (blue, dashed line) are shown. 



fibers is illustrated by Fig. [T] for both types of disorder. 

The sliding of a fiber means that its equilibrium length 
gets suddenly increased, which has the consequence that 
the fiber can only keep a lower load ai = E{e—e\y^, where 
tTi denotes the local load of fiber z, ej^ is its sliding thresh- 
olds and e is the macroscopic strain of the system. Under 
a fixed external load cr, the other fibers of the bundle have 
to overtake the load dropped by the one which has just 
been extended. For the load redistribution we assume 
an infinite range of interaction, i.e. global load sharing 
[1, [in ■ However, it does not imply that the load is every- 
where the same in the bundle. At a given macroscopic de- 
formation e the fibers have suffered a different number ki 
of slidings occurring at different threshold values, hence 
the load they keep will be different: a.i — E{e ~ kis^f^) 
for quenched, and Ui = E[e — e^^ — — ... — e^j^') 
for annealed disorder [i = 1 , . . . , A^) , respectively. (A 
fuse model with a slightly similar microscopic dynamics 
was considered in Ref. where fuses burn out due to 
gradual overheating.) 

We start the analysis with the case of quenched disor- 
der and assume that only a single sliding is allowed for 
the fibers kmax = 1- At the macroscopic deformation e 
the fibers with sliding thresholds ej^ > e are still intact 
and have the local load Ee. However, those fibers which 
have already suffered a sliding e\y^ < e keep only the load 
E{e — e\f^) so that the macroscopic constitutive equation 



follows as 

€ 

a{e) ^ Ee[l- P{Ee)\+ j p{Eei)E {e - ei) dei. (1) 

Note that the integral in the second term is performed 
over the entire loading history of the system. In the inte- 
grand the probability density p of the thresholds occurs 
since the load of the extended fibers depends on the pre- 
cise value of the deformation where the sliding initiated. 
Allowing for two sliding events k^ax = 2, at the macro- 
scopic deformation e those fibers slided exactly once 
whose sliding threshold falls in the range e/2 < ej^ < e, 
while the ones with e^^ < e/2 already suffered two slid- 
ings, hence, a{e) can be written as 

cr(e) = Ee[l- P{Ee)]+ [ p{Eei)E [e - ei) dei 

Je/2 

re/2 

+ / p{Eei)E {e - 2ei) dei. (2) 

For any arbitrary value of k^ax the sliding threshold of 
fibers which are intact or suffered a k number of sliding 
events 1 < fc < kmax fall in the following sub-intervals of 

.rain ^raax~\ 
th T'^th J 

e < ej^, intact, 

e ■ e 

< elf, < — sliding k times, where k < kmax 

k + 1 k 

< ej^ < shding kmax times. (3) 
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FIG. 3: Remaining deformation Sr as a function of the maxi- 
mum deformation e,„ reached before the unloading starts for 
different values of kmax- £r and Em are normalized by the 
average sliding threshold {sth) and the maximum value e™"^, 
respectively. 
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FIG. 4: {Color online) Comparison of the constitutive curves 
of the system with hardening (no failure) and failure after 
kmax sliding events. The arrows represent unloading curves 
for kmax = 2 with a decreasing Young modulus. 



The general form of the macroseopie constitutive equa- 
tion can be obtained by integrating the load kept by the 
above subsets of fibers 

cr(e) ^ Ee[l- PiEe)] 

+ y2 p{Eei)E {e - kei) dei 

k=i Je/(k+l) 

-I- / p{Eei)E [e — ei)dei. (4) 

Note that in Eq. @ the fibers retain their initial stiff- 
ness after suffering kmax sliding events, which is ex- 
pressed by the last term. It can be seen that for small 
deformations £ — > only the intact fibers are relevant 
which result in a macroscopic linear behavior with a 
Young modulus equal to that of the individual fibers. 
For very large deformations e — *■ oo, practically all fibers 
have suffered kmax restructuring events so that only the 
last term of Eq. ^ survives which can be further written 
as 

a{e) Ee - kmaxE j p(ei)ei(iei. (5) 

The integral in the second term converges to the av- 
erage value of the thresholds {sth), hence, we obtain 
ct(£) ^ Ee — kmaxE {eth) ■ The result implies that the 
bundle has an asymptotic linear behavior with the initial 
value of the Young modulus. When unloading the system 
(7{e) — > 0, the fibers simply relax with a linearly elastic 
response since there is not any mechanism in the model 
to decrease the equilibrium length of fibers. It has the 



consequence that as the load is released, the system re- 
laxes along a straight line of slope E equal to the Young 
modulus of fibers and an irreversible remaining deforma- 
tion occurs Er, whose maximum value e™°^ is propor- 
tional to the average slip length (eth) and the number 
of shding events kmax allowed e™"^ = kmax (eth)- Un- 
loading the bundle after reaching a macroscopic defor- 
mation Em the remaining irreversible deformation can 
be obtained as Sr = £m — cr{em)/E. For the purpose of 
explicit calculations we used a WeibuU distribution for 
the sliding thresholds with the cumulative distribution 
P{£th) = 1 — e~'^^"^/^^"' over the range < Eth < +oo, 
where A sets the scale of the threshold values and m de- 
notes the WeibuU exponent. In all the calculations A = 1 
was chosen. The constitutive curve of the model is pre- 
sented in Fig. [2] for the case of m = 1, i.e. exponential 
distribution. It can be observed that increasing kmax the 
asymptotic linear part ofa{e) is preceded by a longer and 
longer plateau regime indicating a plastic response of the 
system. Figure [3] presents the remaining permanent de- 
formation Er of the bundle measured after unloading the 
system from the maximum deformation Em of the loading 
process. One can observe that Er tends to the maximum 
value kmax i^th) as Em increases. 

In realistic situations it is typical that after a certain 
number of sliding events the relaxed length of the sub- 
units of the system cannot be further extended, which 
then results in breaking. Such cases can be captured 
by the model with the assumption that after a fiber has 
slided kmax times it breaks irreversibly, i.e. its stiffness 
is set to zero. The constitutive equation of the sliding- 
breaking system can be obtained analytically by skipping 
the last term of Eq. ^ which has the consequence that 
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for large deformations e — > oo the stress tends to zero 
a{e) — > indicating that the bundle gradually looses its 
load bearing capacity. Figure [H presents a comparison 
of constitutive curves Eq. ^ with hardening and fail- 
ure after kmax slidings. It is important to emphasize 
that due to the breaking of fibers, the unloading Young 
modulus of the system is a decreasing function of the 
maximum deformation e„i achieved before the load is re- 
leased so that the remaining deformation takes the 
former = em-criem)/E[l-P{Eem/kmax)]- InFig.Hthe 
arrows indicate unloading curves where the decreasing 
unloading modulus can also be observed. Under stress 
controlled loading the constitutive curve a{e) can only be 
realized up to the maximum, where macroscopic failure 
occurs. The value of the maximum defines the fracture 
strength ctc ~ cr(ec) of the sliding-breaking bundle, where 
the critical deformation Sc can be determined numeri- 
cally from the equation 1 = £c[^k=T~^i'^/^)Pi-^^c/k) — 
(kmax - 'i-)/kmaxP{Eec/ Knax)] ■ Numerical analysis of the 
above equation shows that both Sc and Cc are increasing 
functions of kmax converging to finite values. 

In the case of annealed disorder, the fibers get a new 
threshold value from the same probability distribution 
after each sliding event (see also Fig. [H6)). The new 
threshold may represent that the physical properties of 
the new conformation attained by the fiber are differ- 
ent from the previous ones including also the effect of 
thermal noise. The macroscopic load a on the system 
can be obtained by summing the load above subsets of 
fibers with different sliding numbers k = 0, 1, ... , kmax 
as a{e) = CTo(e) + cri{s) + ■ ■ ■ + crfe„,„^(£)- Due to the in- 
dependence of subsequent thresholds the terms of the 
above constitutive equation can be cast in the forms 



ai(e)- J p{ei)d{ei)E{e^ei){l- Pie-ei)), (6) 

€ S — €l 

(72 (e) = J J deid£2p{ei)p{£2) (7) 

^th ^th 

X [1 - P(e - ei - £2)] E{e - ei - £2). 



Finally, for k = kmax we get 



e e-ei ek„ 



e 

xE{e- ^ ei), 



Y[ deip{ei) 



1=1 



(9) 



1=1 



Figure [5] shows also examples of the constitutive curve 
a{e) of the fiber bundle with annealed disorder. We 
obtain qualitatively the same behavior with the same 
asymptote as with quenched disorder, however, the plas- 
tic plateau becomes significantly longer. The result 
shows that the extension of the plastic regime is deter- 
mined by e™*^^, however, the precise functional form of 
a{e) is sensitive to the type of disorder. 

We presented a micro-mechanical model for systems 
which can extend their length in a series of stick-slip 
events when subjected to an external load. The model 
is an extension of fiber bundle models in such a way 
that single fibers have stick-slip rheology characterized 
by quenched or annealed threshold values of the local 
strain. We showed analytically that varying its parame- 
ters the model provides a broad spectrum of constitutive 
behaviors: for a large number of sliding events a plas- 
tic regime develops which is then followed by hardening 
(no breaking) and softening (fiber breaking). Unloading 
the system a permanent deformation remains which is a 
monotonically increasing function of the maximum defor- 
mation. If the fibers break after kmax sliding events, the 
unloading modulus goes to zero with increasing defor- 
mation. The constitutive curves provided by the model 
have qualitative agreement with the measured response 
of stick-slip systems [l|, [Hj, however, for a quantitative 
comparison further improvement is needed: allowing for 
an increase or decrease of the Young modulus of single 
fibers after local sliding events would make the model 
more realistic. Assembling bundles of fibers to form a 
chain (serial coupling of fiber bundles) is also a promising 
direction to reproduce the mechanical response of biolog- 
ical tissues. Work in this direction is in progress. 

F. Kun acknowledges the Bolyai Janos research fellow- 
ship of the Hungarian Academy of Sciences. 



The general case of arbitrary k < kmax reads as 



E e—ei e— ei— £2 £k 



1=1 



Y[deip{ei) 

1=1 

k 

E{e-Y,ei), (8) 
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